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Abstract 

Using first-principles calculations, we systematically study the adsorption behavior of a single 
molecular H2O on the Be(0001) surface. We find that the favored molecular adsorption site is the 
top site with an adsorption energy of about 0.3 eV, together with the detailed electronic structure 
analysis, suggesting a weak binding strength of the H2O/Be(0001) surface. The adsorption inter- 
action is mainly contributed by the overlapping between the s and p z states of the top-layer Be 
atom and the molecular orbitals l&i and 3ai of H2O. The activation energy for H2O diffusion on 
the surface is about 0.3 eV. Meanwhile, our study indicates that no dissociation state exists for 
the H 2 O/Be(0001) surface. 
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Water adsorption on metal surfaces has gained a lot of interest in a variety of phenomena 



such as heterogeneous catalysis and corrosion of materials 



, l2j . As a result these systems 



have been intensively investigated by various experimental and theoretical techniques. With 



respect to the substrates, most of the recent studies have 



ocused on the transition metal 



surfaces, such as Cu(100) , Fe(100) [5|,l6|, andPd(lOO) [7|,|8|. Although ambiguities exist, 
by early density-functional theory (DFT) calculations, theoretically a flat-lying configuration 
on the top site of transition metal surfaces has been established by some sophisticated studies 
111 ], arguing that the molecular orbitals (MOs) of water, mainly l&i, dominate the water- 



surface interaction, by coupling with atomic d orbital of the transition metal surfaces. In 
comparison with the vast studies on the adsorption properties of H 2 on transition metal 



12] have 



surfaces, however, for simple metals which are lack of d states, very few researches 
been reported up to now. This long-term overlook should be altered by considering the fact 
that there is now an increasing practical demand for a thorough understanding of the water 
structures on some specific simples metals. 

Motivated by this observation, in the present paper we use first-principles calculations 
to investigate the adsorption properties of H2O on the Be(0001) surface. The reason why 
we choose beryllium as the prototype for simple metals is as follows: (i) beryllium has vast 
technological applications due to its high melting point and low weight. One of its important 



usage is as a plasma facing material in experimental nuclear fusion reactors 



13|]. As the first 



wall, the beryllium can adsorb the residual gases in the plasma vessel, improving the vacuum 
cleanliness [14]]. Water is the main residual gas in the ultrahigh vacuum (UHV) vessels of 
the fusion reactors, so it is highly meaningful to study the adsorption of water molecule at 
beryllium surfaces; (ii) besides and prominently, different from the bulk, beryllium surfaces 
have a large directional s and p electronic distributions around the Fermi energy and thus 
may display unique interaction with the adsorbed H 2 molecules. 

Through analysis of projected density of states (PDOS) and charge density difference, 
we obtain the adsorption properties of H2O on the Be(0001) surface. It is found that the 
H2O molecule prefers to adsorb on the surface top site in a fairly flat-lying configuration 
with a weak binding strength. The adsorption interaction is mainly characterized by the 
overlapping between the s and p z states of the top-layer Be atom and the MOs lbi and 
3ai of the water molecule. We notice that different from transition metal surfaces, the MO 
3ai plays an important role in the interaction between H 2 and Be(0001). The diffusion 
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energetics of H 2 across the Be(0001) surface is also studied, which gives an activation 
barrier of ~0.3 eV. 

Our calculations are performed within DFT using the Vienna ab-initio simulation pack- 
age (VASP) [if]]. The PW91 [isj] generalized gradient approximation and the projector- 
augmented wave potential are employed to describe the exchange-correlation energy 
and the electron-ion interaction, respectively. The cutoff energy for the plane wave expan- 
sion is set to 400 eV. The Be(0001) surface is modeled by a slab composing of five atomic 
layers and a vacuum region of 20 A. A 3 x 3 supercell, in which each monolayer contains 
nine Be atoms, is adopted in the study of the HoO adsorption. The water is placed on one 
side of the slab only and a dipole correction [18| is applied to compensate for the induced 
dipole moment. During our calculations, the bottom two atomic layers of the Be surface 
are fixed, and other Be atoms as well as the H2O molecule are free to relax until the forces 
on the ions are less than 0.02 eV/A. Integration over the Brillouin zone is done using the 



Monkhorst-Pack scheme 



191 ] with 7x7x1 grid points. And a Fermi broadening 20j of 



0.1 eV is chosen to smear the occupation of the bands around the Fermi energy (Ep) by a 
finite- T Fermi function and extrapolating to T = K. 

The calculation of the energy barriers for the water diffusion processes is performed 
using the nudged elastic band (NEB) method 21], which is a method for calculating the 
minimum energy path between two known minimum energy sites, by introducing a number 
of "images" along the diffusion path. Then the energy barrier is determined by relaxing the 
atomic positions of each image in the direction perpendicular to the path connecting the 
images, until the force converges. In present work, the diffusion path is modeled using seven 
images, two of which include the minimum energy sites as initial and final positions, with 
five linearly interpolated, intermediate images between them. 

The structural and energetic parameters of the free water molecule are calculated within 
a box with the same size of the adsorbed systems. The optimized geometry for free H 2 
gives a bond length of 0.97 A and a bond angle of 104.4°, consistent with the experimental 
values of 0.96 A and 104.4° 22]. The calculated lattice constant of bulk Be (a, c) are 2.26 



A and 3.58 A, respectively, in good agreement with the experimental measurements of 2.29 
A and 3.588 A [23J. 

As depicted in Fig. Q^a), there are four high-symmetry sites on the Be(0001) surface, 
respectively the top, bridge (bri), hep and fee hollow sites. H 2 adsorption at all four 
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FIG. 1: (Color online) (a) The structure of the p(3 x 3) surface cell of Be (0001), and three on- 
surface adsorption sites. The red balls denote the initial positions of O atoms in the adsorption 
picture, (b)-(d) Top view (upper panels) and side view (lower panels) of the optimized structures 
of three most stable adsorption states of H2O/Be(0001) surface, i.e., top-rc, top-yl and top-f/2, 
respectively. Gray, red and white balls denote Be, O, and H atoms, respectively. 

high symmetry adsorption sites on Be(0001) surface is considered. The O atom of water is 
initially placed on the precise high-symmetry sites with various orientations of water with 
respect to the substrate. We find that there exist locally stable adsorption states on the top 
site of Be(0001), where the H2O molecules lie fairly flat on the surface, labeled by employing 
the notations top-x, yl and y2, respectively. The structural and energetic details of the 
molecular states are illustrated and summarized in Figs. []Jb)-(d) and Tabled respectively. 
In Fig. W(Jo) we define two angles (ft and 6. (p represents the azimuthal angle of H2O with 
respect to the surface, and 9 represents the tilt angle between the H 2 molecular dipole 
plane and the surface. The adsorption energy of the system is calculated as follows: 

E a d = -E^O/BelTlOOl) _ En 2 o — -E'Be(OOOl), (1) 

where Eu 2 o, -EBe(oooi)> an d -EH 2 o/Be(oooi) are the total energies of the H 2 molecule, the clean 
Be surface, and the adsorption system respectively. According to this definition, a negative 
value of E a <i indicates that the adsorption is exothermic (stable) with respect to a free H2O 
molecule and a positive value indicates endothermic (unstable) reaction. 

From Table [J we can clearly see that at these stable adsorption sites, the work function 
3.75 eV is much smaller than the clean Be(0001) surface (5.41 eV), implying an observable 
charge redistribution between the adsorbate water and the surface Be atoms. Take the 



4 



TABLE I: Calculated structural parameters, adsorption energy for a water molecule on Be(0001) 
surface. E a (eV) represents the adsorption energy. <3? (eV) represents the work function, zq (A) 
represents the vertical height of the O atom from the surface. do-H (A) represents the bond length 
between the O and H atoms. (°) represents the tilt angle between the H2O molecular dipole 
plane and the surface. oh-O-h (°) represents the H-O-H bond angle. 



Site 


E a $ z 


do-H 


9 


«H-0-H 


top-x 


-0.300 3.75 1.79 


0.98 


21.6 


106.3 


top-yl 


-0.303 3.75 1.79 


0.98 


20.7 


106.4 


top-y2 


-0.300 3.75 1.79 


0.98 


20.6 


106.4 



adsorption site top-yl for example, the O-H bond length 0.98 A is almost identical to 0.97 
A of free H2O molecule, but the H-O-H bond angle 106.4° is larger than that of free H2O. 
The tilt angle is 20.7°, which is in line with the value of about 20° of the H 2 O/Al(100) surface 



12J, but differing from that (~10°) of the H 2 0/transition metal surfaces, where the top site 



is also the most stable state. The adsorption energy —0.303 eV indicates a weak molecule- 
surface interaction, and the molecular state may diffuse easily with a low barrier, passing 
through the bridge sites as a transition state, which will be discussed in the latter part of 
this paper. More interestingly, we find that for a H2O molecule to reach the adsorption 
state, there does not exist any energy barrier, which means that H 2 can be adsorbed on 
the Be(0001) surface spontaneously. This is quite different from the adsorption of other 
molecules on the Be(0001) surface such as 2 24 1 and CO |25| , which need to overcome 
energy barriers. 

As presented in Table HJ the energetic differences among these three adsorption states on 
top site are very tiny. For further illustration, we investigate the azimuthal orientation of 
adsorbed H 2 0, which is shown in Fig. [2j It can be seen that the favored orientation for H 2 
at the top site is top-f/1. We can attribute this to the symmetry of the adsorption structure. 
The Be atom at the second layer (on the hep site) of the substrate might also be partially 
responsible for this tiny azimuthal orientation preference. Nevertheless, it is apparent that 
the azimuthal rotation is facile, therefore the adsorption of H 2 on the top site is insensitive 
to the azimuthal orientation. In addition, a considerably less stable adsorption state, where 
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FIG. 2: (Color online) Calculated adsorption energy of H2O as a function of the azimuthal angle 
at the top site. The insets show the structures adopted in the calculations. 

the water molecule is adsorbed on the top site in an upright configuration with the O atom 
down, is found by our calculations, the adsorption energy of this state is about 0.19 eV, 
which is also insensitive to the azimuthal orientation. In the following discussion, we will 
focus our attention on the the most stable flat-lying adsorption sites on the top site. 

In order to gain more insights into the precise nature of the chemisorbed molecular state, 
the electronic PDOS of the H2O molecule and the topmost Be layer are calculated. As a 
typical example, here we plot in Fig. |3]the PDOS for the most stable adsorption configu- 
ration of top-yl. For comparison, the PDOS of the free H2O molecule and clean Be(0001) 
surface are also shown in Fig. |3j Three-dimensional (3D) electron density difference Ap(r), 
which is obtained by subtracting the electron densities of noninteracting component systems, 
p Be (° 001 )( r ) + p H2 °(r), from the density p(r) of the H 2 O/Be(0001) surface, while retaining 
the atomic positions of the component systems at the same location as in H2O/Be(0001), 
is also shown in the inset of Fig. [3]^b). MOs 2a± and 16 2 of water (not shown here) are far 
below the Fermi energy and thus remain intact in water-metal interaction. Here we consider 
only three MOs of water below the Fermi energy, i.e., lb 2 , 3ai and lbi. One can clearly see 
from Fig. [3](a) that after adsorption, these three MOs shift down in energy 3.8 eV, 4.8 eV 
and 4.9 eV, respectively. Compared to free water molecule, the distortion of the adsorbed 
water, i.e., the effect of water geometry, is pretty small (see Tabled]). Hence, these shifts of 
energy can be attributed to the interactions with the substrate. In addition to energy, it is 
noticeable that the MOs of adsorbed water are broaden apparently for lb\ and 3ai, which 
are known to have an oxygen lone-pair character perpendicular to the molecular plane, and 
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a mixture of partial lone-pair character parallel to the molecular plane and partial O-H 
bonding character, respectively Furthermore, despite little change in the PDOS near 
the surface Fermi energy, the water adsorption introduces new peaks for both s and p z states 
of the surface Be atom, aligning in energy with lb\ and 3a!. Especially for 3ai, more elec- 
tronic states of the Be atom appear nearby, compared with lb\. This is quite in accordance 
with our above-mentioned result that the tilt angle of H2O on the Be(0001) surface is larger 
than that on the transition metal surfaces, where only the lb\ dominates the water-surface 
interaction by coupling with atomic d orbital of the surfaces. Moreover, it may also give a 
reason why there exists a less stable adsorption state in an upright configuration on the top 
site. 

All these features indicate a covalent coupling between the adsorbed water and the sub- 
strate. This result is substantiated by the 3D electron density difference. We can see that 
there exists a large charge accumulation between the adsorbate and substrate. As depicted 
by the PDOS, however, no apparent charge transfer between water and surface Be atom can 
be observed. In addition, we find that at the second layer of the surface, there exists a small 
charge depletion right at the fee site, which possibly confirms our earlier speculation that 
the state top-yl is most stable at the top site. 

Given a H2O molecule at a stable adsorption site, it is interesting to see how it can 
diffuse across the Be(0001) surface. Here, therefore, we calculate the diffusion paths and 
energetic barriers of water on Be(0001) surface between neighboring adsorption sites along 
the top-x, top-yl and top-y2 channels respectively, which are schematically shown in the 
insets of Fig. HI The adsorption energies as a function of the lateral displacement of O atom 
are shown in Fig. HJ We can see that the activation energies for H 2 diffusion along the 
three paths are 0.289, 0.291 and 0.286 eV, respectively. The transition state of diffusion is 
located on the bridge site, which is not a stable adsorption site. The similarity of the energy 
barriers between these paths implies that it does not have a strong preference for a particular 
direction in which to diffuse, consistent with our previous analysis illustrated in Fig. |2j By 
assuming the attempt frequency 10 13 of the adsorbate, the energy barrier corresponds to 
diffusion temperature of about 113 K, suggesting that the diffusion can occur under room 
temperature on the H 2 O/Be(0001) surface. 

Finally let us discuss the possibility of dissociation of H 2 on the Be(0001) surface. Ex- 
perimentally, it was found that beryllium can be oxidated by the water vapor [26] . Basically, 
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FIG. 3: (Color online) The PDOS for the H2O molecule and the top-layer Be atom bonded to 
H2O, at the stable adsorption site top-yl: (a) free H2O and the clean Be(0001) surface, (b) the 
molecular state of H2O/Be(0001) surface, where the inset shows the 3D electron density difference, 
with the isosurface value set at ±0.025 e/A 3 . The area filled with yellow color represents molecular 
orbital of H2O. The Fermi energy is set to zero. 

since the binding strength of H2O/Be(0001) system is weak, hence one anticipates that it is 
not prone to dissociate for H 2 on the ideal Be(0001) surface under low temperature. This 
is confirmed by our first-principles static as well as molecular dynamics simulations. We 
find that the hollow hep and fee sites are stable both for H and OH species on the Be(0001) 
surface. When we put the H and OH species simultaneously at two nearest-neighbor hollow 
sites, however, it turns out that they spontaneously combine together into H 2 molecule 
finally without any barrier encountered. This is the same for the next-nearest-neighbor con- 
dition. The MD simulations are performed using the Verlet algorithm with a time step of 1 
fs within the micro canonical ensemble. The O atom of the H 2 molecule is initially set to be 
4 A away from the metal surface. The substrate atoms are initially at rest, while the initial 
kinetic energy of H 2 is set to be 0.6 eV, with the initial velocity towards the substrate. In 



the same conditions it has been found that O2 molecule would dissociate finally [27|. Here 
for H2O, nevertheless, it rapidly reaches a stable state that is similar to the adsorption state 
found above. Consequently, it proves that no dissociation state exists for the H 2 O/Be(0001) 
surface under room temperature. We speculate that the experimentally observed surface 
oxidation of beryllium by water vapor is closely related to the surface roughness or defects. 
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FIG. 4: Diffusion of H2O on the top site of Be(0001) surface as a function of the lateral displacement 
of O atom from its original site top- re (upper panel), top-yl (middle panel), and top-y2 (lower 
panel), respectively. 

In conclusion, we have systematically studied the adsorption behavior of H2O on the 
Be(0001) surface by using first-principles DFT method. It has been found that the water 
molecule prefers to adsorb on the top site by lying fairly flat on the surface, insensitive to 
the azimuthal orientation and with a weak binding strength. It has been shown that the 
MO 3ai, as well as lbi, plays an important role in H^O-Be interaction by overlapping with 
the s and p z states of the underlying beryllium atom. The diffusion energetics of H 2 on 
the top site across the Be(0001) surface has been calculated to display an energy barrier 
of about 0.3 eV. Moreover, no dissociation states of water molecule are found for the ideal 
Be(0001) surface. 
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